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Abstract. Let XW = {X^ ;t > 0}, /i > 0, be the n-dimensional hyperbolic Brownian 
motion with drift, that is a diffusion on the real hyperbolic space H" having the Laplace- 
Beltrami operator with drift as its generator. We prove the reflection principle for X^) 
f*^) | which enables us to study the process X' 1 *' killed when exiting the hyperbolic half-space, 

£N| ■ that is the set D = {x € H" : x\ > 0}. We provide formulae, uniform estimates and describe 

asymptotic behavior of the Green function and the Poisson kernel of D for the process X^\ 
Q ■ Finally, we derive formula for the A- Poisson kernel of the set D. 

r— v 1. Introduction 

PLh ■ In recent years we have seen a considerable growth of interest in stochastic processes on hy- 

perbolic spaces. In a series of papers [I], [2J, [I], [5], [13], [E], [20] many important aspects of 
the hyperbolic Brownian motion were investigated, such as formulae and properties of tran- 
sition density functions, harmonic measures and Green functions. The hyperbolic Brownian 
motion is a canonical diffusion on hyperbolic spaces. It has the half of the Laplace-Beltrami 
operator as its generator. Consequently, it is a hyperbolic analogue of classical Brownian 
motion on Euclidean space. From the other side it is closely related to geometric Brownian 
motion and Bessel processes, see [7] and [6]. Thus, it has some important applications in the 
■ risk theory in financial mathematics, see [9] and [19] . 

The present paper is natural continuation of the research started in [I], where the integral 
representation of the Poisson kernel for the interior of the horocycle, along with the resulting 
analysis of the asymptotic behavior was given. We consider the hyperbolic Brownian motion 
with drift /x > 0, in the half-space model H n = {x E R n : x n > 0} of n-dimensional 

real hyperbolic space, i.e. a diffusion with a generator |A M , where 



(N 



k=l "I* 



Note that for /i = {n — l)/2 the process is a classical hyperbolic Brownian motion 

on H n . We describe the Green function Gd(x,i/) and the density function of the harmonic 
measure (the Poisson kernel Pjj(x,y)) of the hyperbolic half-space, i.e. the set D = {x e 
HI" : X\ > 0}. We provide formulae for those objects as well as the uniform estimates for the 
full range of parameters x and y. Note that the operator A M is a strongly elliptic operator 
on every bounded subset of H n (bounded in hyperbolic metric). Then, the general result 
of [S] implies comparability of the hyperbolic and Euclidean potential theory (comparability 
of corresponding Poisson kernels and Green functions). Since D is unbounded, this kind of 
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result is no longer true and deriving explicit formulae for Pp and Gd seems to be the only 
way to analyze the asymptotic behavior of those objects. It is also worth mentioning that 
the uniform estimates presented in Theorem [3] and Theorem [5] together with those obtained 
in [6] are the only results of this kind in the case of unbounded subsets of hyperbolic spaces. 

The paper is organized as follows. In Preliminaries we provide basic information about 
Bessel processes and hyperbolic spaces. In the part related to Bessel processes we follow 
the exposition given in [T5] and [IB], where we refer the Reader for more details and deeper 
insight into the subject (see also [17]). In Section [31 after describing the structure of the 
hyperbolic Brownian motion with drift we provide formulae for the transition density 

function, potentials and global Poisson kernel for the process In Proposition H] we 

prove the reflection principle which gives the formulae for the transition density function for 
the process Xp \ i.e. the process X^ killed at the boundary of the hyperbolic half-space 
D. In Sections H] and [5] we introduce the formulae for the Green function and the Poisson 
kernel in terms of the modified Legendre functions and examine the asymptotic properties of 
those objects as well as its uniform estimates. All the properties of the Legendre functions 
used in those sections are given in the Appendix. We end the paper with Section [BJ where 
an application of the Girsanov's Theorem together with the results of Section lead to the 
formulae for the A-Poisson kernel of the set D and the solution of the appropriate Dirichlet 
problem on D for the Laplace-Beltrami operator with drift. 



This section contains some preliminary material related to Bessel processes and geometric 
Brownian motion as well as real hyperbolic spaces which are indispensable in the sequel. 

2.1. Bessel process and geometric Brownian motion. We denote by J?M = {R^\t > 
0} the Bessel process with index v G R starting from Rq = x > 0. In the case — 1 < v < 0, 
it is when the point is non-singular, we impose killing condition at 0. Then, for every 
v < 0, the transition density function is given by (see [3] p. 134) 



Here I„(z) is the modified Bessel function. 

Let us denote by B = {B t ; t > 0} the one-dimensional Brownian motion starting from 
and by = {B^ — Bt + ut;t > 0} the Brownian motion with a constant drift v 6 R. 
Let be the geometric Brownian motion with drift starting from x > 0, i.e. 



Bessel process and geometric Brownian motion are related one to the other by the Lamperti 
relation, which gives 



2. Preliminaries 




x, w > 0. 



(2.1) 



Y {u) = {xexp{Bi u) ),t > 0}. 




(2.2) 



where the integral functional Ax (t) is defined by 



A { :\t)= f {Y^\s)) 2 ds = x* [ 



Jo Jo 



exp (2B S + 2vs) ds. 



(2.3) 
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Whenever p is strictly positive the limit of when t goes to oo exists a.s. and its 

density function h^} is given by (see [9]) 

r 2[i -xl/2u 

h ^ {u) = roT)2M Wm) ' (2 ' 4) 

2.2. Real hyperbolic space. For every n = 1, 2, 3, . . . we consider a half-space model of 
n-dimensional real hyperbolic space, i.e. the set HP = {x = (xi, . . . , x n ) G R n : x n , > 0} with 
Riemannian metric 

ds = —• 

The hyperbolic distance dn"^, between two points x and y from HP is described by 

x — y\ 

cosh(du"{x,y)) = 1 + — . (2.5) 

^XnVn 

The group of isometries of HP is generated by reflexions with respect to hyperplanes {x G 
R" : Xi = const.} and inversions with respect to the spheres perpendicular to the hyperplane 
P = {x e K n : x n = 0}. Thus, the set 

D = {x G H n : X! > 0} (2.6) 

can be considered as hyperbolic analogue of the half-space. Indeed, the set D and the interior 
of the complement of D can be transformed one to the other using the isometry HP 3 x = 
(a?!, x 2 , ■ ■ ■ , x n ) — > (—Xi,X2, ■ ■ ■ , x n ) G HP. We will denote by d±D = {x G HP : x± = 0} 
the boundary of £) in H n and by d 2 D = {x G R n : x n = 0, X\ > 0} the second part of the 
boundary of D considered as a subset of R n . Note that d 2 D does not belong to HP. We will 
write dD for d\D U d 2 D. The hyperbolic volume element is given by 

dV n = —dxt...dx n , (2.7) 

X n 

where dx\...dx n stands for the Lebesgue measure on R n . Finally, the Laplace-Beltrami 
operator associated with the Riemannian metric is 

A H"=^E^2 ~{n-2)x n - w -. (2.8) 

k=1 UX k UX n 

This operator is the only second order elliptic differential operator on HP, annihilating con- 
stant functions, which is invariant under isometries of the space. 

2.3. Notation. For every x = (x±, . . . , x n _i, x n ) G R n we introduce the following notation 

x = (xi, . . . , x n -i) G R n_1 , x = (-xi, x 2 , . . . , x n ) G R n . 

Consequently, for every subset A C R n we denote A = {x G R n : x G A}. Moreover, for 
every x, y G HP we put 

p = d m »(x,y), p = d m ™(x,y). 
Note that using ( 12.5ft we easily obtain the following relation between p and p 

cosh p = cos p + — — , 1,1/6 1" (2.9) 

Xnyn 
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Finally, we write / rs g, x G A whenever c~ 1 f(x) < g(x) < cf(x) for every 16 A 

3. Hyperbolic Brownian motion with drift and reflection principle 
3.1. Definition and structure. For every \i > 0, we define the following operator 

n d 2 d 

k=i k n 

Note that for \x = the operator is the Laplace-Beltrami operator (12. 8p . 



The hyperbolic Brownian motion (HBM) with drift is a diffusion = {X t , £ > 0} on 



HP starting from = x G H n with a generator |A^. For /i = we obtain the standard 



HBM on IF. 

The structure of X^' can be described as follows. If (3 = {(/3i(t), . . . ,(3 n -i(t)),t > 0} is 
(n — l)-dimensional Brownian motion starting from (x±, . . . , x n -i) G R™ -1 independent from 
a geometric Brownian motion y(~W starting from x n > then 

= GSiC^W),---,^!^^)),^) (3-i) 

Moreover, following the Lamperti relation one can also show that 

* t (Ai) = (A(i4t M ^0),-»,A.-i(4: M) (0),i2 ( " M) (4; /i) (0)) , (3-2) 

where = {Rf\t > t} is a Bessel motion with index —fj, starting from x n > 

independent from Brownian motion (3. 

3.2. Transition density function, potentials and global Poisson kernel. In this sec- 
tion we compute basic characteristics of the process 

IW suc h as a transition density func- 
tion, its potential kernel and so-called global Poisson kernel, i.e. the density function of 
lim^ooX^. 

Proposition 1. Transition density function (t, x, y) of X^ with respect to the canonical 
volume element dV n is 

W ir(2it)"IVI'-> \ 2 )h (coshu + coshp)'"* 1 " 2 
where v = and p = G?e« (x, y) . 

Proof. The density function f^\ of (A { ^\t), Y^^) was computed in [19] . We have 

fSM = ^^e-^i-exp ("^^) ^/u(t), (3.4) 

where 

r (t) = e^e^) sinh(6) sin d6. (3.5) 

Moreover, the Laplace transform of the function 8 r is (see |18| ) 

oo 

e- x %(t)dt = /^(r), (3.6) 
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where I u (z) is a modified Bessel function. Relation (13. ip and the fact that Brownian motion 
(3 and geometric Brownian motion Y^~^ are independent allow us to write 

/DO 1 / 1 n — 1 \ 

The factor appearing in front of the integral follows from formula ( 12.7ft . Consequently, 
using the formulae ( 13.4)) . ( 12. 5 p and substituting m = x n y n / s we obtain 

p<">(f,s,y)=P) e —— s^e— h "9 s (t)ds. (3.7) 

\VnJ W Jo 

The formula ( 13.51) and Fubini theorem lead to 

poo i>oo (n 2 -b 2 )/(2t) / irh\ f°° 

/ s^e- scoshp 6Jt)ds = / — — — — sinh(6) sin I — / jVe^H^^sfft 
Jo Jo [2iiHyl 2 V t J Jo 

and the substitution s = — t—^ — rr in the inner integral gives 

cosh p+cosh o ° ° 



fi-v p -fi 2 t/2 



e-^l 2 f°° e (- 2 - fc2 )/ 2 * sinh b sm(nb/t) 



y n J 7r(27r)"/2ti/2 J Q ( cos h 6 + cosh p) 



(n+l)/2 



^ v e~^ 2 /n+l\ f°° e^ 2 - 62 )/ 2t sinh6sin(7r6/t) „ 

do. 



^7 \ 2 ;y (cosh6 + coshp) (n+1)/2 

□ 

An easy consequence of Proposition [1] is the following relation between transition density 
functions for the processes with different indices. 

Corollary 1. For every x, y G HP we have 
where v = 



pM(t, x, y) = ( ) exp ( t ) J (t, x, y), (3.1 



Another consequence of the formula ( 13. 3p is the representation of the potential kernel of 
the process which is defined by 



oo 



VM (x, y) = / (t, x, y)dt, x,y eW' 



in terms of modified Legendre functions. 
Proposition 2. For every x, y G HP we have 



(t \ 2e- { ~ v - 1 ' 2 ^ Q U i 7 /o(coshp) 



(2ir)»/2 sinh"" 1 /' 



where v = an d P = <^h™ (x, y) ■ 
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Proof. Representation (13. 7p of p^(t,x,y) in terms of the function 9 r together with the 
formula ( 13. 6p for its Laplace transform imply 

poo / I pco poo 

/ pM(t,x,y)dt=l^) -V/ s ^ e - scoshp e-^ 2 9Jt)dt ds 
Jo \VnJ (^r Jo Jo 



The integral formula (see [12] 6.622 p. 702) 



r s ,-i e - S co Sh p / {s)ds = [i e - {v - w « Q^/^p) 

Jo V vr sinh" 1/2 p 



ends the proof. 



(3.10) 



□ 



Finally, the representation of the Legendre functions in terms of elementary function 
provided in Appendix gives 

Corollary 2. For odd n > 3 we have 

n — 3 

y{p) i n (^Y~ U ^t±A V (coshp-sinhp)^ fc 
1 ,V) \y n ) (27T)" f^r(fi + l + k) (BiDhp)V2+* ' 

where v = p = d^n(x, y) and a{m, k) are the coefficients of the Bessel polynomials 

(m + k)\ 

aim, k) = - — — — — — . 

v ' ' (m-k)\k\2 k 

Proof. Follows immediately from (17. 7p and (13. 9p . □ 

We denote by P^(x, y) the global Poisson kernel, i.e. the density function of the random 
variable X$ = lim^oo X^\ where x = Xff^ and y G R n l x {0}. The existence of this limit 
follows from the representation (13. ip and the fact that Y} ^ tends to zero a.s. whenever 
H > 0. 

Proposition 3. For x G H n , y G R n_1 x {0} we have 

r(pi + v) xf 

\x - y\ 

where v = 



pMfr y) - f3ir 



Proof. Using the representation (13. ip . the fact that = a.s and the formula (12. 4p we 

obtain 

1 / i~ -i 2 ' 

1 I \x — y\ 



pM(x,y) = / — == exp - l - — ^- h^\u)du 



x 2 / f°° ( \x-y\ 2 + x 2 



exp . 

r(/x)2^(y2¥) n ~ 1 Jo \ 2u 
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Making the substitution the last integral we arrive at 

pM(x, y) = ^— — I" e~ v v^ v ~ x dv - + v) 



r(^(\x-y\ 2 + xlY +U Jo r0*)7r" \x-y\ w »- 

□ 

3.3. Reflection principle. We denote by r D = inf{£ > : X^' ^ D} the first exit time 
of the process X^*) from the hyperbolic half-space D (see (P for its definition) and 
stands for the process killed when exiting D. Its transition density function is given by 
the Hunt formula 

p%\t 1 x,y)=pM(t,x,y)-&[t<T D wM(t-T D ,Xto 1 y)], t > 0, i^gD. (3.12) 

We introduce the following reflexion principle for HBM with drift. Recall that using the 
notation introduced in the Preliminaries we have X&) = (— X[^\ X^f \ . . . ,Xn). 



Proposition 4 (Reflection principle). Let us define 



XF>, r D <t. 



then {Z[%>» i {*<">} 



t>o- 



Proof. Fix m — 1, 2, . . ., = t < *i < ... < t m < t m +i = oo and 4 . . . , A m G B(W). 
Then, for every k — 0, 1, . . . , m we have 

P x (t k <r D < t k+l , 4 M) 64..., e A m ) = 
= P x (h <r D < t k+1 , zg> 64-.., 4 M) e A, z £i e A k+i, • • • > ^? e A«) 
= P x (4 <t~d< tk+i, X t f 6 4,-.., -X'tf G 4 X t ^ +l 6 4+i> . . . , 6 A m ) 
Applying strong Markov property we obtain that the above-given expression is equal to 

W(t k <r D < t k+1 , 64-, g A k , E x 'o (x%l_ Tn 6 4^,..., xtl TD G 4)) 

Note that the first coordinate of Xt£ is equal to zero a.s. and for every t > we have 
p(o,x 2 ,...,x n )( X M e A) = p(o,*2,...,«n)(x t M G A). Thus, we can change 44.., 4 into 
A k+ i, . . . , A m and usage of the strong Markov property leads to 

P x (t k <r D < t k+u zg> 64,.-, Z% G 4) 

= P x (t k <r D < t k+l , xg> 6 4 ••• ; Xj£ 6 4) 

Summing up the above-given equality with respect to k we obtain that the finite dimensional 
distributions of the processes 1^ and are the same and it ends the proof. □ 

The consequence of the reflection principle, as in the real case, is the following formula for 
transition density function of the killed process. 

Theorem 1. For every x,y G D and t > we have 

p ( g\t,x,y)=pM(t,x,y)-pM(t,x,y). (3.14) 
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Proof. Whenever A G B(D) and t > we can write 



P x {Z < t ) e A) = P x (t < r D ; e A) + P x (t > t d] g A) 



= P x (t < r D ; X^ eA) + P x (X { / l) e A) 

= P x (t<r D ; Xf EA) + P x {X ( f 1 G A). 
Consequently, applying Proposition HI we obtain 

p ( g\t,x,y)dV n (y) = [ p^(t,x,y)dV n (y)- f p^(t,x,y)dV n (y) 



p%\t,x,y)dV n (y)- / p%\t,x,y)dV n (y)- 

A J A 



□ 



4. Green function 
The Green function of the set D is defined as usual by 



poo 

G%\x,y)= / p { £\t,x,y), x,y G D. 
Jo 



As an immediate consequence of Theorem[T]we get the following representation for (x, y) 
in terms of the potential kernel V^'{x,y). 

Theorem 2. For every x,y G D we have 

G ( g (x, y) = VV> (x,y)- (x,y). (4.1) 

In the Propositions OlS we examine the asymptotic behavior of G£ (x, y) at the boundary 
of the set D and whenever x — y tends to 0. 

Proposition 5. Let z G d±D, x G D. Then 

l^f 1 DK ' (27r)^/ 2 x n z n \z n J (smhd Mn (x,z)y+^ K > 

Proof. Using the representations (14.11) . (13. 9p . the mean value theorem and the formula (17. 5p 
we obtain 

G%\x,y) = vM(x,y)-VM(x,y) 

QlZi/l( cosh p) _ Q^j/2( cosh p) 

4e -(v-i/2)in xm / Xn y- V (sinhp)^- 1 ^ ~ (sinhp)"- 1 ^ 



(2ir) n / 2 x n y n \y n J cosh p- cosh p 

v+l/2 
H-l/1 



4e— '*x iyi fx n y- v QTvM 



(2tt)-/ 2 x n y n \y n J (^~l) n l 2 
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where coshp < 9 < cosh p. Using the fact that the function z — > e~ mn ^ 2 Q l/ ^ 1 y 2 2 (z) is 
decreasing on (1, oo) we get 

4e-"/ 2 C!g (coshp) (x n \ v -"x n y n m 4e~™ ' 2 ggjg (coshp) 

(2tt)™/ 2 (sinhp)"/ 2 - Vy„7 arij/i D ~ (2tt)"/2 ( s i n hp)"/ 2 ' 1 j 

By (j2.9p we get lim^^g^ coshp = coshp and consequently, by continuity of the Legendre 
function, we obtain (14.21). □ 



Proposition 6. Let z G d 2 D, x G D. Then 

lim ylT n G^ (x,y) = ^ "\ { - - , %^ ) . (4.4) 



2/S-D 



X — Z X — 2 



Proof. Using ( 13. 9 p and (12. 5p we get 

2e -(^-i/2)MT Qrirt (coshp) 

</„. -^n (27r)n/2 y n sinh ^i/2 p 

^ m^h ( 1 y +u ^ coch ^ ^ :i //2(coshp) \ 

The asymptotic formula ( 17.4ft gives 

r(p + z/) x 2 / 



\imy l - n V^\x,y) - 
l^D ^V(ii+l)\ x - z \^ +2v 

Changing y to y in above given equality and using (14. ip we prove (14. 4p . □ 

Proposition 7. For x G D, (y\, ...,y n -i) G R"" 1 we have 

lim y^G%\x, y) = AV{ ^ V ^% l^ iyi . (4.5) 

Proof. Multiplying (14.31) by y£ +u+1 we obtain 

Ay^Ql% (coshp) ^^+1^+2 4^+^ (coshp) 

(2tt)«/2 ( s inhp) n /2 - xxyxe-™/ 2 D {X ' V} ~ (2vr)™/ 2 (sinh p)"/ 2 ' K '> 

By formula for hyperbolic distance (12. 5p we have 



lim -jg-= lim ( - + 1 = 2x n . (4.7) 



= lim I — + ± - y 

y»->-oo COSh p fe->oc \ y n 2x n ?/„ 

Rewriting the right-hand side of (14. 6p in the form 



2^ - 1 



4 / v V + " +1 ( (coshp)' 
' ^ » 'Wip)^ 1 



(27r) n / 2 V coshp / l v F; (sinhp)^ +1 / 

and applying ( I7.4p and ( 14. 7p we get that it tends to 



4r(p + t/ + r 

7T I T(p+ 1) 



-rwri/2 n+v+X 
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as y n — > oo. By (12. 9p we obtain lim^-^ coshp/coshp = 1 and the limit of the left-hand 
side of (14. 6p is exactly the same as above and consequently (14. 5 P holds. □ 



Proposition 8. For fixed x G D we have 

1 71-2 



lim \x 



>J\ 



' D 



x,y) 



X 



n-2 



2W 2 ' 



Proof. Whenever x G D (i.e. X\ ^ 0) we have x ^ x and \im y ^ tX V^\x^y) = V^(x,x) 
which is finite. Consequently 



lim 

y^-x \x 



VW(x,y) 



|2-n 



lim 

ii-KB \x — y 

2g-invr/2 



2-n 



\im(2x n y n 



(2vr) n / 2 r-*r 

— inn/2 

lim(x n y n 



n — 2 
2 



a—v 



\x - y\ 



n-2 
2 



Qm-i/2( cos1i p) 
(sinh p)^- 1 /2 



n-2 
2 



— 1 (cosh p 

Vn, 



1/2 (cOSh p) 

J (sinhp)^- 1 ^ 



.r 



n-2 



2vr n/2 



r 



□ 



where the last equality follows from (1T.3H . 

We end this section by the following uniform estimates of the Green function Gp(x,y). 
Theorem 3. There exists constant c = c(n, p) > such that 



r 

Vn 



)1 — U 



\%nyn 



■ I/-1/2 



\x — y\ \x — y\ n ~ (coshp)^ +1 / 2 
where p = dj$n(x,y), v = V£ ^- 

Proof. Using formulas (14. ip . ( 13. 9 p and applying Lemma [3] we get 



x,y G D, 



x 

y 

x 

y 

2^+3/2 



p—v 2g-(^-l/2)i7r 

(2tt)™/ 2 



sinh^ 1 / 2 



sinh 1 '" 1 / 2 



P 



a—v 



cosh p — cosh p 



2/n 



cosh p — 1 (cosh p — l) i/_1 / 2 (cosh p)^ +1 / 2 
x x y\ (x n y n ) u - 1/2 



\x — y\ \x — y\ n (cosh p)^ +1 / 2 



□ 



5. POISSON KERNEL 



The random variable X<f£ is supported on d\D whenever Td < oo and if td — oo then 
G d<iD. We denote by Pp(x,y) the Poisson kernel of D, i.e. the density function 
of random variable Xr^ with respect to the Lebesgue measure (on d\D and ^-D). In the 
following theorem we provide explicit formula for Pjy(x,y) in terms of Legendre functions 
and global Poisson kernel P^(x,y). 
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Theorem 4. For every x G D and fi > we have 



11 



2xi y 



[l—U—l 



(2tt)"/ 2 x ^ v 



x 



[ r(/i)7r- V|x-y 



2/J 



1 2^+2^ 



2 M \ 



(5.1; 



where v = p = dm™(x,y). 



Proof. Since takes values in <9 2 -D if and only if Td = oo, for A G BipzD) we get 

P x (Ig 6A)=P J (Xg> E A) =P X (t d = oo; Xj' G A) 
= P x (JW G A) - P* (re < oo ; G A) . 



Applying (13. lip and the reflection principle given in Proposition H] we arrive at 



x 



2fi 



2\^+ v 



(\ x - y\ ) 

2fi 



dy 



X 



2/' 



T(fJ,)ir v 



x: 



\x-y 



2^+2u 



r(/x)7r" 7i (|x-y| 2 ) 

2/L4 



dy 



\x — y\ 



To deal with the case X-j-n G 5i-D, (i.e. td < oo) recall the representation (13. 2p . We denote 
by ri = inf {s > : /?i(s) < 0} the first exit time of the one dimensional Brownian motion 
Pi from the positive half-line. Since the function t i — > A x ~ (t) is continuous, increasing and 
^4x^^(0) = a.s., we have Ai„ (td) = T\ whenever tjj < oo. Thus, for every C G B(diD), 
we can write 



P* (iW G C) = P x (td < oo, I T W G C) 
= p(t 1 <^(oo),(/3 
= P(r 1 <4r^(oo),(/3 n , J R^ ) ) GC). 



p(-M) 



G C 



Moreover, we have Ax n (°°) = W > where T V w = inf{t > : i?- 



(-m) 



0} is the first 



hitting time of zero for Bessel process R( ^\ i.e. the lift-time of R( Since (f3 Tl ,R 



n 



C C diD implies R^ T1 ^ > and consequently we can remove the condition n < A Xn (oo) 
and write 

P* (Xg G C) = P* ((/3 n) ^: M) ) e C) = P X ((^(r!),...,^-!^),^^)) G C) 
The density function of n is given by 



»7*i( u ) 



2,71 



u > 0, 
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and since T\ and the process (/3 2 , . . . , P n -i, R^~^) are independent we have 



n-l 



P%\x,y) = [ — — ^ziexp ( ^- V(x fc - y k f ) g ( u ^\y n )r] Xl (u)du 
Jo (27rn) n/2 \ 2u tZ J 

=#5* (£)T exp (- £ ? cosh ") (^) """ +s,/w 

Here the last equality follows from (12. ip and (12. 5ft . Substituting u = x n y n /s and then using 
( KTUh lead to 

™ „./j,—u—1 roc 

^ ) (x,y) = 2(27r)TO 1 /2 „ 1/2 ^ + , I exp ( - S coshp)/,( S )^ 1 )/^ 

r^ 1 /2 __ /2 Q^ 2 1/2 (coshp) 



g l 2/n / — mri/2 M 1/ 

( 27r )n/2-l/2 ^ VtT sinh"/ 2 



(2tt)«/2 sinlr^p ' 

□ 

The explicit formula for Pp{x,y) presented above leads to the uniform estimates of the 
Poisson kernel P^\x,y). Note that the constant appearing in the estimates given below 
depends only on n and the /i > 0. 

Theorem 5. There exists c = c(n, fi) > such that for every x G D, y G dD we have 



P%\x,y) 



XlXnyl* 1 „ 

V\ = o, 



i \ n (\ i 2 i V 
— 2/| +x n y n ) 



\x 



2/1 (1^ — 2/| J 
where v = and p = d^n (x, y) . 

Proof. It is an easy consequence of formulas ( 15. ip . ( 17. lip and ( 17. 9p . 



□ 



6. A - Poisson kernel 

In this Section we will use the notation to denote expectation with respect to the 
distribution of the hyperbolic Brownian motion with drift /i > starting from x G HP and 
consequently, we will drop the superscript (/i) in the notation of the process whenever it 
appears under E^ } . 

For any bounded domain U CC H n we introduce the following probabilistic definition of 
the A- Poisson kernel (A > 0) for HBM with drift /i > 

P(}'\x, y) = EM [e- Xrv -X TU G dy] , xeU,ye dU . (6.1) 

where tjj is the first exit time of the process from the set U . In the following theorem 
we introduce the connection between A-Poisson kernels and Poisson kernels for HBMs with 
different drift parameters. 
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Theorem 6. For any bounded domain U CC H" and fi,r] > we have 

E(m) 



13 



p-T) 



E?) [X r , G dy] . 



(6.2) 



Proof. The last coordinate of HBM with drift is of the form xjfi (t) = x n exp ( Wt — rjt) , where 
W t = B n (t) + (77 — /j)i for one-dimensional Brownian motion B n . According to Girsanov's 
theorem, for every T > 0, process {W t , < t < T} is a standard Brownian motion with 
respect to measure defined by 



dQ 



T 



dP 



exp (/j - rj)B n {T) 



(77 - 



T 



and consequently, the process jx^, < t < Tj is HBM with drift with respect to measure 
Qt and its generator is \A n . Hence, for any A C U we obtain 

E?) [tu < T; X TU G A] = E«* fa < T; X^ G A] 

U < T; exp ( (fj, - V )B n {t v ) - ^ ~^ Tu) ]X Tu eA 



k < T; (exp (S n fa) - fiTu))^ e ~W~^ ; X T[/ G A 



7t/ 



< T; (X n (n,))^e-3^-^; G A 



Since U is bounded we have tjj < 00 a.s. and by monotone convergence theorem, taking 
limit as T goes to infinity we obtain 



E?) [X^ g A] = <-"E s M [^fa))^-^ 2 -^^ G A 
which is equivalent to (16.21) . 



□ 



In the case of unbounded sets U (such as hyperbolic half-space D) the first exit time tjj 
can be infinite with positive probability and then the right-hand side of ( 16. ip vanishes. From 
the other side, taking the sequence of bounded sets U m CC HP with smooth boundary such 
that U m /* D and using the formula (16. 2p we can easily see that the A-Poisson kernel P£ ,X 

multiplied by where 77 = \/2\ + tends to Pp . This observation leads us to consider 
the following Dirichlet problem on D. 

For given / G Cb(dD) and A > we are looking for u G C 2 (D) satisfying 



-A,.u (x 



\u(x] 



, x G D, 



(6.3) 



such that Xn 2X+ ^ M u(:r) is bounded and 



lim x^/ 2X+ ^ ^u(x) = f(z), z G 3D. 



(6.4) 



Consequently, we extend the previously given definition of A-Poisson kernel for bounded set 
U into the case of unbounded hyperbolic half-space D. 
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Definition 1. For any A > 0, the A-Poisson kernel of the set D for HBM with drift fj, > 
is a function P£ x : D x 3D — > [0, oo) such that for given / G Cb{3D) 

u ( x ) = / f(y)PD X ( x ^y) d y 

JdD 

satisfies conditions ( 16. 3 p and ( 16. 4p . 

Justification for the correctness of the definition as well as the formula for Pj£ x are given 
in the following 

Theorem 7. The function P£' x described in Definition^ is unique and it is given by 



D 

(2tt)"/ 2 ^+ 1 sinh^p 
T(a + u) ( x« +fl a% +li 



y e 3 2 D, 



(6.5) 



\x — y\ \x — y\ 



T(a)-K u 

where v = r\ = y? + 2A. 

Proof . For a given / G Cb(3D), let « be a C 2 function on .D satisfying (16. 3p and ( 16. 4p 
and consider the sequence of bounded sets U m CC H n with smooth boundary such that 
U m /* D. For every m, there exists a unique solution w m G C 2 (U m ) to the above-given 
Dirichlet problem with D replaced by U m and / = u G Cb(3U m ), see [TT|. Moreover, it is 
given by 

%(i)=Sf ) [e-^-u^J]. 

Since tt is also a solution, using ( 16.21) we get 

u(x) = EW [e" A ^ M (X Tl7 J] = <-"E(") [(X^ m )^« (X T[7 J] . 

By dominated convergence theorem we get 

u{x) = xJT'EW [/ {X TD )\ . (6.6) 

It gives the uniqueness of the solution of the above-given Dirichlet problem and consequently, 
it proves correctness of the Definition [TJ If we now define u by (16. 6p . then obviously h(x) = 
x 1 l~ fl u(x) = [f (X TD )] is bounded and the condition (16. 4p holds by continuity of sample 
paths of the process. Since A v h(x) = 0, we have 

iA,u(x) = l -\ {x^>h(x)) 

= x^-A.Hx) + v)<- v+1 ^r(x) + ^^^Kx) 

= <~ V \ (a m - 2(/x - rf)x n ^j h(x) + (/i - ri)x^ +1 -^-(x) + \x^h{x) 

= ^M(x) + AxJ-Mx) = A«(x). 
Finally, the formula (16. 5p follows from ( 16. 6 p and Theorem HJ □ 
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7. APPENDIX 

Since many formulas appearing in the paper involve the Legendre functions Qf(z) for 
convenience of the Reader we collect here some basic information and properties of those 
functions. Here we follow the exposition given in [TU] (see also [12])- In the second part of the 
section we introduce several lemmas providing some detailed properties regarding differential 
formulas and uniform estimates. 

The Legendre functions Ql{z) and P£(z) are solutions of the Legendre's differential equa- 
tion ([12], P-958) 

(1 - z 2 )u" - 2zu' + (a(a + 1) - b{l - z 2 )' 1 ^ = 0, b + 1 > a > --. 

For a > 0,b G R the asymptotic behavior of the function Q b (z) is described by (see [10], 
3.9.2 (10), (21)) 

lim (z- l) a/2 Q a b (z) = e am 2 a/2 ^T(a). (7.1) 



limz «« w = v^x^;,: 1 ) . ( r.2) 



2 b+1 r(6 + 3/2 

Using the above-given properties of Q%(z) we easily obtain 



„ +b+1 Q a b (z) _V^e am T(b + a + l) 



lim z a+b+1 b = v L (7 A) 

z-XL Z (^_l)a/2 2 b +!r(6 + 3/2) ' 1 1 



whenever < a < b. Finally, we introduce the following 
Lemma 1. Let z > 1. Then for a,b G R we have 

d Qt(z) Ql + \z) 



dz(z 2 -l) a / 2 (z2_l)(o+l)/2> 

Proof. We apply equality (see [10], 3.8 (1) p. 160) 

(6 + a )(b - a - l)Ql\z) = Q a b +1 (z) + 2a^== Q a b (z) 

y z z — 1 

to the formula formula (see |10| . 3.8 (9) p. 161) 



and we get 



Hence 



(7.5) 



|qg(*)= (fl + ^ + V H*)-^QgW 

l z Ql(z) = QV\z) + -£L_ Qt{z) . (7.6) 

d Q a {z) _ (z 2 - l)^ 2 j- z Qt(z) - azQt(z)(z 2 - l) a / 2 ~ l 
dz~(z 2 - \) a l 2 ~ (z 2 - l) a 

_ (z 2 - 1)§-Q a b (z) - azQ a b (z) _ Q a + i {z) 

= ( z 2 _ iy/2+l = ( 2 2 _ 1 )(a+l)/2' 
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where the last equality follows from (17. 6p . □ 
Lemma 2. For every m — 0, 1, 2, 3... and b G R we have 

^ (Z)=Z( " 1} V2 r l +m+ 2j^ r(6 + 3/2 + fc) ^g^ W^' (7 - 7) 

where a(m, k) are the coefficients of the Bessel polynomials 

(m + fc)! 

a(m '* )= (m-fc)!fc!2*- 

Proof. To prove the result we apply mathematical induction with respect to m. For m = 
and 6 G R the equality (17.71) is just the formula 3.6.1 (12) in |10j . Assume thesis for some 
m > 1 and any 6 G R. Denoting 

(z - y^M )" 
(z 2 - 1)« 

we have following recursion 

d 



/?(*) = ,1 — n V € R,g > 0, 



^/ 9 P W = -(P + 2g)/ 9 P +1/2 (.) - 2 g yj£(*). 
It comes from bellow-given calculation 



d (z - V^iy- 1 M z (z - Vz^iy 

Tz^ z) =p ^ - 1). (1 - 7P = T ) ~~ qz 

= ~ P {Z 2 _ 1)g+ l/2 " 2 ^ - (2 2_ 1)ff+ l + 



(z* - 1) 



9+1 



According to (17. 5p and inductive assumption we get 

"2 ? (-i) m+1 gr +1+1/2 (^) /2 *(-i) m 9 g™ +1/2 ( 



r (6 + m + |) ^3T( m+1 )/ 2+1 /4 V 7T T (6 + m + |) ^ v /^_TT m+1/2 

/ (m + fc)! , v 

~ 9 * (m-k)\k\2 k T(b + 3/2 + ky^+V*' 



m (m+fc)! 

■ E fffp^ ((» + m + | + a)^ W + (m + i + W 

fc=0 * 

■ g (m-fc)!fc!2*r (& +f + fc) (6 + m + 2 + 2A;)/ f 

v^ 1 (m + fc - 1)! 

(m _ fc + i)!^ _ !)!2fc-ir(6 + | + k) [m + K)I f+i+r Z 
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m + | 3 6+1 (2m)!(2m + 2) 3 



r(6+f) v 2' 2 +iV ' m!2 m r(6 + § + m) y 2 +1+^+1)2 



fe=i 



(m+fc)! (m+fc)! \ 

(m-fc)!fc!2 fc , , 3 . (m-fc+l)!(fc-l)!2 fe -i \ ,,b+±+fc 

_^ (6 + m + _ + 2i) + __^_ 



We simplify the expression in last parenthesis as follows 

(m+fc)! (m+fc)! 
(m-fc)!fc!2* , A 3 (m-fc+l)!(fe-l)l2*-i 

T(b+l + k) [b + m+ 2 +2k)+ r(6 + i + Jfc) 



(m+fc)! 

(m_fc+1)!fc!2fc ' (6 + m + ^ + 2k){m - k + 1) + 2k{b +\ + k) 



T(b + l + k) V V " ' "" ' 2 ' ™"' v '" " ' ~' ' ""^ 1 2 

(m+fc)! 

jpp|( 6 + m + ?)( m + , + 1 ) 

(m+fc+1)! 
(m-fc+l)!fc!2 fc 3. 



r((, + § + *) (i, + m + 2- 

Consequentely 

^ i(-l) m+1 QT +1+1/2 (z) 
71 T (6 + m + |) /^j( m+1 )/ 2+1 /4 ~ 



& + + m + 3) M (z) _ £5^ + 2) + 3 - + m ) f b+ ^ +m ( 

r(6 + |) ( + + 2 j/ f + lW T(b+\ + m) {b+ 2 + m)/ f + i + ^ + D 2 ( 



2; 

(m+fc+l)! 



E(m-fc+l)!fc!2 fe 3 b+\+k , . 
^ - o + m+- r„ , fe m 

fc=l v & ' 

m+1 (m+fc+1)! 
-_(J + m + !)V (m-fc+l)!fc!2* b+I+fc 

m+l (m+fc+1)! _ ^-x 6+ i +/ r, 

/ 7 , 3.v^ (m-fc+l)!fc!2 fe Z-VZ^-l +2 + 

= — (o + mH — > ir r e — • 

2 S r (&+§ + *0 _i)f+i+| 

Multiplying polar expressions by (z 2 — iy m / 2 +3/4 /i;_[_!_Jl^) we g e ^ f ormu l a j n thesis for 



2 

m+l. □ 



Lemma 3. For 1 < x < y we have 



(x 2 - l) a / 2 (y 2 - l) a / 2 



c(a,b) y — X 



y-l(x- l) a x b+1 ' 



{11 



Proof. For x > 1 function is continuous and positive. Using equalities ( 17.3ft and ( 17.41) 

we obtain 

-q«r £ 1 C7 

( x 2_ 1 )a/2 (x-l)^ 1 1 J 
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Let assume now that (y — 1) < 2{x — 1). Mean value theorem and Lemma [T] gives us 

QtM Q a b (y) 



c - ai7T (x 2 -!)"/ 2 fa 2 -l)°/ 2 = e _( a+ l)i 7r _ 



Qt +1 (0 



y-X (^2_ 1 )(a+l)/2' 

where £ £ y]- The right-hand side function is decreasing, hence 



-(a+l)i7r 



(£2 _ l)(a+l)/2 



< e 



-(a+l)i7r 



(a; 2 - l)(«+i)/2 - (x-1) 



a+l x b+l 



< 



2c 



y-l(x- l) a x b+1 
In case of (y — 1) > 2(x — 1) it holds that 



Than 



^2 _ 1 )a/2 ( y 2 _ X)a/2 



y — x 
2- > 1. 

y-i ~ 



inni 2c 

< 



(7.10) 



Ql( x ) JZ3J c 



[x 2 - \y/ 2 

y-x 



(x - l) a x b+1 



y — 1 (x — l) a x 



Ctrr-b+l ' 



To prove the another inequality we estimate 



Otto 



Ql(y) 



_( x 2_ 1 )a/2 ( y 2_ 1 )a/2 



y-i 

y-x 



-i x— 1 

2/-1 



(x - l) a x' 



a 3+1 



(x-l) a x 



x — 1 



a~6+l 



6+1 



(y - i) a y 



aM Qt(y) 



jm i 

> 



i - 



x-l 

y-i 



"1 


fx — 








4)1 


c 


- c \y- 





The last expression tends to - as 2— r — )■ 0. There exist e > such that for ^4 < e 

r c y-1 y-1 



( y 2 - iyi 2 



x-l 



Ql(x) Qt(y) 



_(x 2 -l) a / 2 (y 2 -\y/ 2 \ \y-x 



x - l) a x b+1 ) > — 



2c 



what is equivalent to 



(x 2 - l) a / 2 (y 2 - l) fl / 2 



> 



1 y — x 



2cy-l(x- l) a x b+1 ' 



For 2— r > e we use mean value theorem theorem again 



« (x 2 -i)°/ 2 (y 2 -i)°/ 2 = p -( a+ i)fa (0 > c -( a +i)fa <?6 + (y) 



y - x 



(£2 _ l)(a+l)/2 



(y2 _ l)(a+l)/2 



1 

> - 



> 



c(y- l) a +V +1 ~~ c {y-l)(x- l) a x b+1 ' 



□ 
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Lemma 4. For 0<x<y,a>0we have 

1 1 c y — x 



x a y a yx c 



Proof. Assume that y > 2x. It follows that 1 > > ~. Hence 

y — x 1 



1 1 




1 






< 




< 2 


x a y a 




x a 




1 1 


> 


1 




x a y a 


x a 






y x a ' 



11 2 a - 1 y - x 1 
> £ 



If y < 2x we use mean value theorem 



_L 1_ 



y -x 

where z G Then we get 



a — TTi 



1 1 2 

< < 



^a+i x aJr ^ yx a 

1 111 

> > 



(7.11; 



ya-\-l 2 a yx a 

□ 
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